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Abstract 

Spin-dependent deconfined interaction in the QQ system is derived 
from the field correlators known from lattice and analytic calculations. 

As a result hyperfine splitting is found numerically for charmonium, 
bottomonium and strangeonium in the range T c < T < 2 T c . Spin- 
orbit interaction due to magnetic correlators (the Thomas term) is 
able to produce numerous QQ bound states with accumulation point 
at M = rriQ + tuq. Possible influence of these effects on the thermo¬ 
dynamics of quark-gluon plasma is discussed. 

In honor of A.Di Giacomo on his seventieth birthday. 


1 Introduction 

The Spin-Dependent Interaction (SDI) in the QQ system is being studied for 
the last three decades ra-ca for heavy qurkonia in the confinement phase. 
The first estimates of SDI essentially exploited perturbative expansion 
and nonperturbative SDI appeared in the form of the Thomas term JJ. In 
the general approach Pi using 1/Mg expansion for heavy quarks, four 
spin-dependent static potentials E(r),i = 1,2,3,4 are expressed through 
connected field correlators. 
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In this approach perturbative expansion of correlators reproduces per¬ 
turbative part of Vi(r) while the nonperturbative part can be obtained only 
from lattice simulations, done by several groups umn 

In the Field-Correlator Method (FCM) |9j (see dS for a review) the 
SDI potentials have been obtained using 1 jM q expansion in JTTj, TSj, in the 
form containing only Gaussian correlators D(x), D i(x), and allowing for 
immediate check of Grornes relation [3J. The nonperturbative component 
of Vi{r) is calculated directly through D(x), D\(x) and the Thomas term 
naturally appears as asymptotics at large distances. Extensive and thorough 
lattice calculations of correlators D(x),Di(x) done by the Pisa group HI] 
have made it possible to determine numerically all 4 potentials V). 

The Erst comparison of resulting spin splittings of charmonium and bot- 
tomonium levels was done in J4J and has shown a good agreement of lattice- 
based potentials, with experiment. At the same time spin-orbit splittings 
were found to be sensitive to the behaviour of held correlators D(x), D\(x) 
at small x. 

In [To] a detailed analysis was done of SD potentials, obtained from held 
correlators in m , as compared to the definitions of Eichten-Feinberg-Gromes 
(EFG) 011. In the latter case SD potentials V) are not shown to satisfy 
Grornes relation [0, and analysis of ra has demonstrated the difficulties 
arising in this respect for the lattice-dehned V) BUl- 

it was understood later M, that the same set of SDI potentials I/)(r) 
can be obtained without the 1 /M q expansion, using cluster expansion of 
gluonic helds for quarks of arbitrary mass m q and in this case instead of 
quark masses rn q , there appear in V)(r) the so-called einbein masses p q , with 
effective average values (p q ) playing the role of constituent masses calculated 
through m q and string tension a. 

This approach works well in the confinement phase for meson m and 
baryon (HI states of light and heavy quarks, when SDI can be considered 
as corrections, see e.g. the analysis [TH] for light scalar mesons, and |Q2IJ for 
heavy quarkonia. 

The inclusion of chiral dynamics (where SDI is crucial) needs modification 
in the formalism, and instead of SDI potentials one considers in this case the 
effective Lagrangians, as in [20], and the huge mass splitting between p and 
7T is obtained as the chiral nonlinear amplification of the standard hyperhne 
SDI. 

So far so good for the zero temperature QCD. For T > 0 there ap¬ 
pears a hrst change in SDI namely one must distinguish colorelectric and 
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colormagnetic correlators, contributing to V)(r), since instead of two Gaus¬ 
sian correlators D(x), Di(x) one has for T > 0 four independent correlators: 
D E (x), D H (x), Df(x), Df{x) and it was agrued in El El 123] that only one 
of them vanishes for T > T c . These correlators have been found numerically 
on the lattice [23] , and the original expectation was confirmed, that the cor¬ 
relators do not change significantly for T <C M^, where is the dilaton 
mass of the order of 0 ++ glueball mass, ~ 1 -y 1.5 GeV. Moreover it was 
shown in [21] that indeed only one of correlators, D E (x ) vanishes exactly at 
T = T C . 

Already at this stage one can show that all four SD1 potentials are ex¬ 
pressed through magnetic correlators, and only one SDI term is due to the 
colorelectric correlator, contributing to spin-orbit force. 

The situation becomes even more interesting for T > T c , where the main 
colorelectric correlator D E (x ), ensuring confinement, disappears but all other 
correlators stay intact. 

In this case all four SDI potentials should change only slightly above T c 
and ensure the spin splitting of levels of the same type as it was below T c . 
Of special interest is the Thomas term, which is doubled in magnitude (due 
to disappearance of colorelectric static force s(r)) and is dominant at large 
distances. It is argued below in the paper, that the spin-orbit Thomas term 
can create infinite number of bound QQ states with the accumulation point 
near 2 Mq. 


2 Spin-dependent interaction in the qq sys¬ 
tem 


For the quark-antiquark Green’s function ( of both light and heavy quarks) 
one can use the Fock-Feynman-Schwinger Representation (FSR) J2S] with 
the kernel containing vacuum fields and quark spin operators in the form 

EU Q2>! 

(W a (x,y)) = (expig J d'K^(z)F flu (z)) = 


where 


= exp^) 


72—1 


(wY 


n\ 


dn(l)... / d7r(n)((F(l)...F(n))), 


dTr^iz) = ds /JjU (z) - ia^dr 


( 1 ) 

( 2 ) 
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Here ds lw is the surface element, while dr is the proper time integration, 
which can be connected to usual time t via dr* = LEI- The spin-spin 
interaction can be obtained from ©, m keeping only quadratic (Gaussian) 
correlators (FF), 



and spin-orbit interaction arises in m from the products (cr pv F pv ds p \F p \) 

imiiEi. 

It is clear, that the resulting SDI will be of matrix form (2 x 2) x (2 x 2) 
( not accounting for Pauli matrices). If one keeps only diagonal terms in 
&nvF pu (as the leading terms for large /i, ~ M) one can write for the SDI the 
representation of the Eichten-Feinerg form 


V. 


(diag) 


SD 


(R) 


(T i L\ a 2 L 2 1 de 2dVi(R) 
l/jf ~ If^^RdR + RdR 


d 2 Li - <?i L 2 1 dV 2 dia 2 V A (R) (3d 1 Rd 2 R - <ji<j 2 R 2 )V 3 

2 fjLiji 2 R dR 12 Hi/j, 2 12 / u 1/ u 2 i? 2 

At this point one should note that the term with in © was obtained 
from the diagonal part of the matrix (m — D)a pu F pu , namely as product 
idkDk ■ &iEi, see (111 Ifi] for details of derivation, while all other potentials 
Vi,i — 1,2, 3,4 are proportional to correlators (Hi$H k &). One can relate 
correlators of colorelectric and colormagnetic helds to D E , Df, D H , D E as 
follow (see pi. m 


g 2 (Ei(x i, (x 2 , r 2 )<f> + ) 


N r 


Sn D E { A, u) + Df + u- 


• dDf 

du 2 


dDf 
+ UiUj -foR 
(5) 


g 2 {H i {x)$H j (y)$+)=N c 


Sij I D h { A, v) + D E + u 


H 


> dP? 

du 2 


UiUj 


dD[ [ 

du 2 


( 6 ) 


BD eh 

g 2 (E i (x)$H j (y)$ + ) = -N c e ijn u n u 4 ^ ) (7) 

As a result one obtains the following connection between SD potentials and 
correlators. 
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1 dVj. 
R dR 


l_dV2 

R dR 


V A 


dv 


r R dX 

'o ~R 


1 - 


V 

R 


D (X, v) 


dv 


rR XdX 


R 2 


D H (X, v) + D E (A, v) + A 2 


H 


dD { 
~d>? 


V, = 




roo ( f)D H \ 

J dv\3D H (R,v)+3D?(R,v) + 2R 2 -^-\ 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 


1 de(R) 
R dR 




D E (X, v) + -Df (A, v) + (A 2 + v 2 ) 


dDf 
dv 2 


( 12 ) 


One can check, that the Gromes relation jT acquires the form 


v;(r)+£'(r)-v'(r) 



[ R dX{D E { A, v) - D h ( A, v)) + l -R{Df{R) - 
Jo 2 

(13) 


For T — 0, when D E = D H , Zlf = D E , the Gromes relations are satisfied 
identically, however for T > 0 electric and magnetic correlators are certainly 
different and Gromes relation is violated, as one conlcl tell beforehand, since 
for T > 0 the Euclidean 0(4) invariance is violated. 

To conclude this section, one comment on nondiagonal terms in ©, which 
contribute to the total Hamiltonian H as m 


D?m 


H = Hofa, fj, 2 ) + V^hh + ^( 75 ) 1 ( 75)2 + ... 


(14) 


where dots imply terms proportional to 11 ( 75)2 and ( 75 ) 112 - It is clear that 
Vgp Tld) has the structure similar to 0 with replacement Vi —■> V- n °" d) and 
contains terms proportional to electric correlators Eq.0. E.g. in [2SJ the 
term yj nond> was found to be 

v} nond) (R) = f dv ( [ 3D E + 3D E + [3v 2 + R 2 )^jPj ( 15 ) 

For heavy quarks (and for light quarks in the states, where SDI is repulsive) 
the extremal values can be found from the minimum of the spinless 
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Hamiltonian H 0 (/j, 1 ,^ 2 ), and in this case SDI gives spin corrections, which 
are not large even for light quarks and are in good agreement with experiment 
- see IH1I2] for heavy and light quarkonia respectively. 

For light quarks in the states with attractive spin-spin interaction the qq 
Hamiltonian should be taken with the full matrix structure as in ED, and 
the stationary values of the quark ’’constituent” masses fif* should be found 
from the minimum of eigenvalues of H , see ESI for more details. In what 
follows we shall consider only diagonal part of SDI, which is essential for 
heavy quarkonia. 

3 Spin interactions in the deconfinement phase 

In the confined phase SDI in terms of held correlators was studied in nu 
U 21 m~m- For T > T c the correlator D E (x ) vanishes, as it was argued 
in HU 123 EH and confirmed in lattice calculations m This fact leads 
to a serious change both in SDI as well as in the spin-independent part 
Ho of the total Hamiltonian H. The latter change was studied in detail in 
recent papers HU 123, where it was shown that in the total static potential 
V(R) = Vn(R) + V\(R), generated by D E (x ) and D E (x) respectively, only 
the term V\ (R) survives for T > T c . This term in contrast to the confining 
potential V E (R ) ~ crR, R —► 00 , saturates at large R and can support 
bound states of cc and bb in some interval of temperatures T c < T < T d , 
while bound states dissociate at T > T d . (T d ~ 2 T c for Df calculated 
analytically in ra, and Td ^ \.\ ( 2!T } C for D^ from lattice calcnlations done in 
m , see for details). All this dynamics is due to colorelectric correlator 

D E and colormagnetic helds were not taken into account. We now consider 
the role of SDI, which is mostly due to colormagnetic helds. 

We start with small distances and assume, according to HU HHl 123 the 
following form of Df , D E , D H 

D em = Df’ H (pert) + D?’ H (nonpert ), D H = D E onpert + D E rV (16) 
where 

16 ry p~Me,. 

Di' H (pert) = - — s -e~ XE ’ HX +0(a 2 s )] Df' H (nonpert, M E , H x > 1) = C e ,h~ - 

Sttx^ x 

(17) 

and ^ 

Dnonpert&Mo > 1) = d H e ~ M ° X , D E £rt = + O(^) (18) 
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Here C'e,h, dn are known constants depending on cr H - the spacial string 
tension, a H = \ J D* onpert (x)d 2 x. 

One can easily check in turn that the small distance behaviour of SDI 
potentials Vi(R) does not change much, since it is defined by the perturbative 
parts of correlators which are little modified when temperature grows above 
T c . The main difference comes at large distances, where according to Eqs. 
mm one finds in the limit of large R, 


Vi(R) = -o„, Vi(R) = ^, 7 h = /_” do £ XdXD^iX, v) ( 19 ) 

| rOO 

e'(R) = - duRDf(R) -> 0, V 3 (R) -»• 0, V 4 (R) -> 0 (20) 

2 .7 — 00 

where in (ED all three potentials are exponentially small at large R. There¬ 
fore the spin-orbit part of V^ ag) has the form (for equal quark pole masses) 

V^(R -* oo) = ^ + 0(e~ MR ). (21) 

The Erst term on the r.h.s. of m was quoted (without derivation) in [2Tj, 
where it was suggested that being dominant at large R (where V\ (R) ex¬ 
ponentially approaches to the constant limit V) (oo)) this term by itself can 
support bound states of heavy (and, possibly, light) quarkonia. Indeed, con¬ 
sidering the Erst term in m as a Coulomb-like potential one arrives for 
heavy quarkonia (where /io ~ m) to the mass spectrum of bound states for 
SL > 0, J = L + 1, 


M n 


2m — 




4 m 3 (L + n + l) 2 ’ 


n = 0,1, 2, 


( 22 ) 


We conclude this section with discussion of hyperfine interaction in the 
deconfined phase. As it is seen in ED, V 4 (R) depends only on D H , Dj 1 and 
does not change across T c according to the lattice data pg. Keeping only 
(the dominant) perturbative part of V 4 , one obtains 


V hs = bm f 17 ^(R). 
9^lh2 

For the hyperfine energy shifts this gives 


A E h f = 


4a a {V'(R)) 

9(^i + jl-2 ) 


-3, S = 0\ 

+ 1 , s=i)- 


(23) 


(24) 
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To calculate (V'(R)), one can use the static potential, computed in |27j, 
m both analytically and on the lattice. A rough estimate can be found 
from comparison of free static energies found on the lattice (see e.g. 0S| and 
refs, therein) both below and above T c (up to T ~ 1.3T C ), which shows a 
close similarity of V'(R ) in both temperature domains. Hence one expects 
for T c < T < 1.3 T c the mass gap between J/ip and rj)c of the same order 
as for T = 0, i.e. 5E ~ 0.1 GeV, which roughly agrees with lattice MEM 
computations [SO) . 

The situation with hyperhne spin splittings in the light qq system is how¬ 
ever different and cannot be treated by the methods given above, since the 
restoration of chiral symmetry at T > T c needs the effective Lagrangian 
technic m mentioned above. The physical reason for that lies in the fact, 
that for light quarks A Ehf Eq. (12-11) . becomes dominant (for = n 2 ~ m) 
and cannot be treated as a perturbation, and one needs to solve nonlinear 
equations for the effective mass operator, given in |20j . This point will be 
treated elsewhere. 

In summary, we have derived the spin-dependent potentials for the qq 
system for any temperature T valid in the situation, when spin splittings can 
be considered perturbatively. 

For 1.3 T c > T > T c spin splittings of charmonium and bottomonium are 
shown to change little compared to zero temperature case. The color mag¬ 
netic confinement produces the Thomas spin-orbit term, which dominates at 
large distances, in absence of colorelectric confinement at T > T c , and can 
possibly support a sequence of bound states. These features demonstrate 
the importance of strong interaction in the quark-gluon plasma, which was 
advocated in [2T]-[21] and supported by explicit calculations in |23 El in 
agreement with lattice data m- 

This work is supported by the Federal Program of the Russian Ministry 
of industry, Science and Technology No.40.052.1.1.1112, and by the grant for 
scientific schools NS-1774. 2003. 


References 

[1] A.De Rujula, H.Georgi and S.L.Glashow, Phys. Rev. D12, 147 (1975); 
T.de Grand, P.L.Jaffe, K.Johnson and J.Kiskis, Phys. Rev. D12, 2060 
(1975). 


[2] T.Appelquist, R.M.Barnett and K.D.Lane, Ann. Rev. Nucl. Part. Sci. 
28, 387 (1978); 

W.Buchmueller, Phys. Lett. 112, 479 (1982). 

[3] E.Eichten, F.Feinberg, Phys. Rev. D23, 2724 (1981). 

[4] D.Gromes, Z.Phys. C26, 401 (1984). 

[5] W.Buchmuller, Y.J.Ng and S.-H.H.Tye, Phys. Rev. D24, 3003 (1981); 
J.Pantaleone, Y. J.Ng and S.-H.H.Tye, Phys. Rev. D33, 777 (1986). 

[6] M.Campostrini, Nucl. Phys.B256, 717 (1985); 

C.Michael and P.E.L.Rakow, Nucl. Phys. B256, 640 (1985); 

P.de Forcrand and J.D.Stack, Phys. Rev. Lett. 55, 1254 (1985); 
M.Campostrini, K.Moriarty and C.Rebbi, Phys. Rev. Lett. 57, 44 
(1986); 

C.Michael, Phys. Rev. Lett. 56, 1219 (1986); 

I.J.Ford, J.Phys. G 15, 1571 (1989); 

A.Huntley and C.Michael, Nucl. Phys. B270, 123 (1986). 

[7] K.D.Born et al. Phys.Lett. B329, 332 (1994); 

G. S.Bali, K.Schilling and a.Wachter, Phys. Rev. D55, 5309 (1997), ibid. 
D56, 2566 (1997). 

[8] M.Koma, Y.Korna, H.Wittig, hep-lat/0510059 

[9] H.G.Dosch, Phys. Lett. B 190, 177 (1987) ; 

H. G.Dosch and Yu.A.Simonov, Phys. Lett. B 205, 339 (1988) ; 

Yu.A.Simonov, Nucl. Phys. B 307, 512 (1988). 

[10] A.Di Giacomo, H.G.Dosch, V.I.Shevchenko, Yu.A.Simonov, Phys. Rept. 
372, 319 (2002); hep-ph/0007223. 

[11] Yu.A.Simonov, Nucl. Phys. B324, 67 (1989). 

[12] M.Schiestl, H.G.Dosch, Phys. Lett. B209, 85 (1988). 

[13] A.Di Giacomo and H.Panagopoulos, Phys. Lett. B285, 133 (1992) ; 
M.D’Elia, A.Di Giacomo, and E.Meggiolaro, Phys. Lett. B408, 315 
(1997) ; 

A.Di Giacomo, E.Meggiolaro and H.Panagopoulos, Nucl. Phys. B483, 


9 




371 (1997) 

M. D’Elia, A.Di Giacomo, and E.Meggiolaro, Phys. Rev. D67, 114504 
(2003); 

G.S.Bali, N.Brambilla, A.Vairo, Phys. Lett. B 42, 265 (1998); 
E.Meggiolaro, Phys. Lett. B451, 414 (1999) . 

[14] A.M.Badalian, V.P.Ynrov, Yad. Fiz. 56, 239 (1993). 

[15] A.M.Badalian, Yu. A.Simonov, Phys. At. Nucl. 59, 2164 (1996). 

[16] Yu.A.Simonov, QCD and Topics in Hadron Physics, Lectures at the 
XVII International School of Physics, Lisbon, 29 September -4 Octo¬ 
ber,1999, hep-ph/9911237 

Yu.A.Simonov, Spin-interactions of light quarks, preprint ITEP 97-89 
(unpublished). 

[17] Yu.A.Simonov, Phys. Rev. D65, 116004 (2002). 

[18] A.M.Badalian, Phys. At. Nucl. 66, 1342 (2003); 

A.M.Badalian, B.L.G.bakker, Phys. Rev. D67, 071901 (2003). 

[19] A.M.Badalian, V.L.Morgunov, B.L.G.Bakker, Phys. At. Nucl. 63, 1635 

( 2000 ). 

[20] Yu.A.Simonov, Phys. Rev. D65, 094018 (2002); 

S.M.Fedorov, Yu.A.Simonov, JETP Lett. 78, 57 (2003). 

[21] Yu.A.Simonov, JETP Lett. 54, 249 (1991). 

[22] Yu.A.Simonov, JETP Lett. 55, 605 (1992). 

[23] Yu.A.Simonov, Phys. At. Nucl. 58, 309 (1995); hep-ph/9311216; 

N. O.Agasian, JETP Lett. 57, 208 (1993). 

[24] M. D’Elia, A. Di Giacomo and E. Meggiolaro, Phys. Rev. D 67, 114504 
(2003); 

A. Di Giacomo, E. Meggiolaro and H. Panagopoulos, Nucl. Phys. B 483, 
371 (1997); 

A. Di Giacomo, E. Meggiolaro and H. Panagopoulos, hep-lat/9603018. 

[25] Yu. A.Simonov, J.A. Tjon, Ann. Phys. 300, 54 (2002); hep-ph/0205165 
Yu. A.Simonov and J.A.Tjon, Ann Phys. 228, 1 (1993). 


10 





[26] Yu.A.Simonov, Phys. At. Nucl. 68, 709 (2005). 

[27] Yu.A.Simonov, Phys. Lett. B619, 293 (2005), hep-ph/0502078. 

[28] A.Di Giacomo, E.Meggiolaro, Yu.A.Simonov, A.I.Veselov, 

hep-ph/0512125 

[29] Yu.A.Simonov, Phys. At. Nucl. 69, N3, (2006), hep-ph/0501182 

[30] P.Petreczky, hep-lat/0409139 


11 






